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Abstract 

Lattice QCD with Wilson quarks and a chirally twisted mass term represents 
a promising alternative regularization of QCD, which does not suffer from 
unphysical fermion zero modes. We show how the correlation functions of the 
renormalized theory are related to the theory with a standard parameterization 
of the mass term. In particular we discuss the conditions under which these 
relations take the same form as obtained from naive continuum considerations. 
We discuss in detail some applications and comment on potential benefits and 
problems of this framework. 



1 Introduction 



Lattice QCD with Wilson quarks Q is widely used to compute hadronic ob- 
servables and matrix elements from first principles. It is a gauge invariant 
regularization with an ultra-local action and an exact global flavour symme- 
try, but all axial symmetries are explicitly broken by the Wilson term. The 
latter fact is usually not considered a fundamental problem, as chiral symmetry 
can be restored by introducing appropriate counterterms. Well-known exam- 
ples are the additive quark mass renormalization, the renormalization of the 
non-singlet axial current and the mixing pattern of the AS = 2 effective weak 
Hamiltonian. The corresponding scale independent renormalization constants 
can be determined both in perturbation theory and non-perturbatively, by 
imposing continuum chiral Ward identities as normalization conditions |2]||. 

As a result chiral symmetry is restored up to cutoff effects, and the prob- 
lem has thus been solved from a field theoretical point of view. However, 
the absence of an exact chiral symmetry has further consequences in practi- 
cal applications, as it implies that the Wilson-Dirac operator is not protected 
against zero modes. This is not a problem in principle, as the functional in- 
tegral over Grassmann variables cannot diverge. After integration over the 
quark fields, a small eigenvalue of the Wilson-Dirac operator appears both in 
the fermionic determinant and in the quark propagators entering the corre- 
lation functions. Fermi statistics then implies that the limit of a vanishing 
eigenvalue is always regular. Despite this fact, numerical simulations with 
some of the standard algorithms may still experience technical problems. In 
particular, one may suspect that accidental zero modes are at the origin of 
long autocorrelation times which have been observed in numerical simulations 
with the hybrid Monte Carlo algorithm 

A conceptual problem arises in the so-called quenched approximation, 
which consists in neglecting the fermionic determinant. The contribution of a 
small eigenvalue to a fermionic correlator is then not balanced by the determi- 
nant, leading to large fluctuations in some of the observables which completely 
compromise the ensemble average ||. Gauge field configurations where this 
happens are called "exceptional" and various recipes of how to deal with them 
have appeared in the literature |]-||] . Strictly speaking, the quenched approx- 
imation with Wilson fermions is ill-defined, as the absence of zero modes is 
only guaranteed for rather heavy quarks. We emphasise that this problem is 
common to all lattice regularizations with Wilson type fermions. However, its 
practical relevance in a given physical situation depends on all the details of 
the chosen lattice action, and on the statistics one would like to achieve in 
numerical simulations. For example, with non-perturbatively 0(a) improved 
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Wilson quarks |],|| and the standard Wilson plaquette action, the problem is 
felt when the quark masses become somewhat lighter than the strange quark's 
mass 0. If the quark mass is further decreased, the frequency of (near-) 
exceptional configurations strongly increases. As the problem becomes even 
more pronounced with increasing lattice volume S, it is clear that the ap- 
proach to the chiral limit with Wilson type quarks is limited by the zero mode 
problem rather than by finite volume effects. 

To solve the aforementioned practical and conceptual problems we pro- 
pose to add a non-standard mass term to the (improved) Wilson quark action. 
The lattice Dirac operator for two quark flavours then reads 

AmQCD = D W + m + iflq'JsT 3 , (1.1) 

where D]y denotes the massless Wilson-Dirac operator, ijiq is the standard 
bare quark mass, and fi q is referred to as twisted mass parameter. It couples 
to a term with a non-trivial flavour structure (r 3 is a Pauli matrix acting 
in flavour space), and protects the Dirac operator against zero modes, inde- 
pendently of the background gauge field. This lattice action has previously 
appeared in a different context [11], and it has already been proposed as a 
regulator for exceptional configurations, implying, however, a limiting proce- 



dure /u q —* at the end [12]. In contrast, our proposal, first presented in [13], 
is to interpret this theory as an alternative regularization of QCD with two 
mass degenerate quarks. We will refer to it as QCD with a chirally twisted 
mass term, or twisted mass QCD (tmQCD) for short. Indeed, in the classical 
continuum limit an axial rotation of the quark and anti-quark fields relates the 
tmQCD action to the standard QCD action. Furthermore the axial rotation 
of the fundamental fields induces a mapping between composite fields. One 
may hence think of this transformation as a change of variables which leaves 
the physical content of the theory unchanged. 

In this paper we want to demonstrate in which sense these classical con- 
siderations can be elevated to a relation between the renormalized correlation 
functions of tmQCD and standard QCD. To this end we first regularize both 
tmQCD and standard QCD using Ginsparg- Wilson quarks (sect. 2). In this 
framework the bare correlation functions can be related by a change of vari- 
ables in the functional integral. Renormalization will be discussed in sect. 3, 
in particular we identify renormalization schemes which preserve the relations 
between the bare correlation functions. Based on universality, it will then be 
clear how to proceed if the regularization does not respect chiral symmetry, 
and we discuss in detail the case of Wilson quarks (sect. 4). We conclude with 
a few remarks concerning current and future work on tmQCD (sect. 5). 
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2 Twisted mass QCD and Ginsparg Wilson quarks 

We start with classical continuum considerations, and then discuss the regu- 
larization with Ginsparg- Wilson fermions. In particular we use a formulation 
which hides the lattice peculiarities as much as possible, so that naive con- 
tinuum relations carry over essentially unchanged to the lattice regularized 
theory. 

2.1 Classical continuum theory 

We consider the continuum limit of the twisted mass QCD actionQ, 

Sj?[il>,$] = j d 4 x^(l/)+m + i W5 T 3 ) if>, (2.1) 

where = + denotes the covariant derivative in a given gauge field 
Gu, and r 3 is the third Pauli matrix acting in flavour space. The axial trans- 
formation 

if/ = exp(ia7 5 r 3 /2)-0, $ = V? exp(ia7 5 r 3 /2), (2.2) 

leaves the form of the action invariant, and merely transforms the mass pa- 
rameters 

ml = mcos(a) + /iqsin(a), (2-3) 
f/ q = — m sin(a) + /i q cos(a). (2-4) 

In particular, the standard action with fj,' = is obtained if the rotation angle 
a satisfies the relation 

tana = ^ q /m. (2-5) 

Chiral symmetry of the massless action leads to the definition of the Noether 
currents, 

^ = V^75yV>, *£ = ^yV, (2-6) 

which are only partially conserved at non-zero mass parameters. More pre- 
cisely, the so-called PCAC and PCVC relations take the form 

d^A* = 2mP a + in q 5 3a S°, (2.7) 
dj^_= -2f, q e 3ab P b , (2.8) 

1 We adhere to the convention that 7-matrices are hermitian, {7^,7^} = 25^, with 
//,!/,... ranging from to 3, and we set 75 = 70717273. 
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where we have defined the pseudo-scalar and scalar densities, 

^ a = ^75y^, S° = ^. (2.9) 

The axial transformation of the quark and anti-quark fields induces a trans- 
formation of the composite fields. For example, the rotated axial and vector 
currents read, 

r a , fcos(a)A? + e 3ab sin(a)K b (a = 1,2), 
\ A n ( a — <j), 

t, t», fco S (a)F; + e 3 <*sin(a)Aj (a = 1,2), 

V >-*TT* = {v„3 (o = 3 >, (2 - U) 

and similarly, the rotated pseudo-scalar and scalar densities are given by 

\P a (a = 1,2), 

P' a = I 1 (2.12) 

[cos(a)P 3 + isin(a)±S° (a = 3), 

S'° = cos(a)S° + 2isin(a)P 3 . (2.13) 

It is easy to verify that the rotated currents and densities satisfy the PC AC 
and PCVC relations (2.7,2.8), with the transformed parameters ml ( |2.3j ) and 



fi' q (2.4). In particular these relations assume their standard form, 



d^A'l = 2m'P' a , dpV'l = 0, (2.14) 



if a is related to the mass parameters as in eq. (|27 

Finally, we note that the tmQCD and standard QCD actions are exactly 



related by the transformation (2.2) and therefore share all the symmetries. 
However, in the chirally twisted basis the symmetry transformations may take 
a somewhat unusual form. For example, a parity transformation is realized by 

ip(x) — ► 7oexp(ia7 5 T 3 )V>(x , -x), (2.15) 
ip(x) — ► ip(x , -x) exp(ia7 5 r 3 )7o, (2.16) 

and similar expressions can be obtained for the isospin and the remaining 
discrete symmetries. It is then straightforward to infer the behaviour of com- 
posite fields under these transformations. 

2.2 Ginsparg- Wilson quarks 

We now replace continuous Euclidean space time by a hyper-cubic lattice with 
spacing a, and we choose some standard lattice action for the gauge fields. The 
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precise choice will not be important in the following, but for definiteness we 
may take Wilson's original plaquette action. As for the quark fields we assume 
that the lattice Dirac operator satisfies the Ginsparg- Wilson relation |14j] , 

D l5 + l5 D = aD l5 D. (2.17) 

This relation arises naturally in the construction of fixed point actions JD]], and 
an explicit solution for D has been given by Neuberger 16]. In the following 
we assume that D is a local operator |0] which satisfies eq. ( |2.17j ) and has 
the conjugation property = 75D75. It then follows that the matrix IS], 

75 = 75(1-^), (2.18) 
is hermitian and unitary. The massless action of a quark doublet 

S F = a 4 ^2$D^, (2.19) 



has a global SU(2) x SU(2) invariance [19] which can be parameterised by the 
transformation 

<// = exp(^r a /2)exp(^7 5 r a /2)^, 

$ = $ exp(i^7 5 r a /2)exp(-^r a /2). (2.20) 

Here Wy^ (a = 1, 2, 3) are real parameters and r a are the Pauli matrices acting 
on the flavour indices. Note that these transformations have the same form 
as in the continuum, except for the appearance of 75. To mask this difference 
we follow refs. 2C|,18] and define left handed fields 

tfL = |(l-W, = ^^(1+75), (2.21) 

and analogously the right handed fields with the complementary projectors. 
The massless action splits into left and right handed parts due to the identity 

$Di/> = faDfa + ^R^R, (2.22) 

and the transformation rules for the chiral fields are exactly as in the contin- 
uum. In particular it is straightforward to construct composite fields which 
transform among each other under a chiral transformation. Examples are the 
isosinglet scalar and the isovector pseudo-scalar densities which may be defined 
through 

S° = i^a + ^ = #-^% (2.23) 
P a = ^ly^r-^ry^l = ^75y(1-5°W- ( 2 - 24 ) 
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One finds 



S' u = S° [i//, ip'] = cos(uj a )S + 2i sm(u A )u a A P a , (2.25) 

where uj a denotes the modulus of the vector (uo\, lo a ,lo\), and u\ = uj\/uj a is 
a unit vector. Similarly, P' a can be expressed as a linear combination of S° and 
P a , and we have thus found a multiplet of bare lattice operators transforming 
under the lattice symmetry in the same way as their continuum counterparts. 
To define the Ginsparg- Wilson regularization of tmQCD we now use the scalar 
and pseudo-scalar densities for the mass terms, i.e. we write 

S F = a 4 [^Dt/} + mS° + 2i/x q P 3 ] . (2.26) 

From the preceding discussion it is then clear that the transformation 

ip' = ex.p(ia%T 3 /2)ip, ip' = ^exp(ia7 5 r 3 /2), (2.27) 



does not change the form of the action, as it corresponds to the special choices 
\ — and u> A = a S 3a in eq. ( ^20|) . Its effect can thus be absorbed in a 



change of the parameters as in eqs. (|2.3|j2.4j ), and the standard QCD action is 
recovered if a satisfies the relation ([2.5]). 

2.3 Functional integral 

We are now prepared to define the functional integral of lattice regularized 
tmQCD with Ginsparg- Wilson quarks. Denoting the sum of the Wilson pla- 
quette and the fermionic action by S, the partition function is defined by 

Z = J D[i)^]D[U]e~ S '. (2.28) 

The physical content of the theory can be extracted from its correlation func- 
tions, i.e. n-point functions of the form 

(0(x 1} x n )) = Z- 1 [ D[?P, ^}D[U] e- S 0( Xl , x n ), (2.29) 



where 0(x\, . . . product of local gauge invariant composite fields 

which are localised at the space time points x\, . . . ,x n . In the following we 
will sometimes indicate the functional dependence upon the quark and anti- 
quark fields in square brackets, i.e. we set O = 0[ip,ip]. 

If one restricts the theory to a finite physical space time volume, e.g. by 
imposing periodic boundary conditions, the functional integral at fixed lattice 
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spacing is a well denned finite dimensional integral, and it is then possible 
to perform a change of variables of the form ( 2.27] ) . As the transformation 



matrices are special unitary, the Jacobian is unity. Due to the form invariance 
of the action the whole effect is a transformation of the parameters in the 
action and of the composite fields, viz. 

{0(x u x n )) M = (0'( X1 , . . . , x n )) {m , K) . (2.30) 

Here the primed field is implicitly defined through 

O'ty'^'] =OhM, (2.31) 



with the transformed quark and anti-quark fields in eq. ( |2~27l ). The subscript 



of the expectation values reminds us of the transformations of the parameters 
in the action according to eqs. ( 2~lj| , 2~4| ). Eq. ( 2.30| ) is an exact identity and 



defines the starting point for the statement of the aforementioned equivalence 
between the renormalized theories. 

In order to prepare the ground for the discussion of the renormalization 
procedure we proceed a little further. First we notice that the combination 

(m') 2 + (^) 2 =m 2 + ^ 5 (2.32) 

is left invariant by the chiral rotation, so that it is convenient to use polar 
mass coordinates, 

m = Mcos(a), / u q = Msin(a), (2.33) 

with radial mass M and the angle a chosen according to eq. (^). Second, 
as the bare theory has an exact chiral symmetry, it is useful to decompose 
operators in multiplets which transform irreducibly under the chiral flavour 
symmetry transformation. We will denote the irreducible components of the 
composite fields by 4>^\(x), where r labels the representation, k identifies the 
multiplet and A labels the members of the multiplet. The general SU(2) x 
SU(2) transformation ( |2.20| ) is then represented by a matrix R^ r \ viz. 

= flS(«v;w A )*SM. (2-34) 
Simple examples are the multiplets of the scalar and pseudo-scalar densities, 
(±S°,iP a ), and of the currents (^l7mT^ l ' V'RTmT'V'r)- 

Without loss of generality we may restrict attention to n-point functions 
of such fields, and the identity ( 2.30| ) then takes the form 

(M,a) 

tl^---€tM) {M ^ al) , (2.35) 
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where we have used polar mass coordinates in the subscript, and the shorthand 
notation 



# (r) (a) = i? (r) (O;0,0,a), (2.36) 



for the rotation induced by the axial U(l) transformation ( [2. 27 ). Note that 



we may keep the angle a' independent of a, and with the choice a' = a the 
r.h.s. of eq. ( |2.35| ) consists of correlation functions with the standard QCD 
parameterisation of the mass term. 

2.4 Ward identities 

Symmetries in Quantum Field Theories are usually expressed in terms of the 
Ward identities which follow from infinitesimal symmetry transformations. 
This has the advantage that identities are obtained between correlation func- 
tions which are defined within the same theory, even in the presence of explicit 
breaking terms. To derive the Ward identities we start by introducing the 
space-time dependent variations of the quark and anti-quark fields, 

SX(u>Mx) = u{xYi{^){x), 6%(u)i>(x) = $(z)-y 5 %u(x), (2.37) 
5^(u)iP(x) =lj(x)^j(x), 5^{u)^{x) = -i>{x)^uj{x), (2.38) 

with some real- valued function u(x). The action on arbitrary composite fields 
is defined by treating the variations like ordinary derivatives. The variation 
of the action ( p. 26 ) then yields 



5l(u)S = -a 4 ^w(x){a^(x)-2mP a (x)-^ q (5 3a 5 (x)},(2.39) 

X 

S$(u)S = -a 4 Y,u(x){d;V«(x) + 2v q e 3ab P b (x)}, (2.40) 

X 

with the pseudoscalar and scalar densities as defined above, and where the 
divergences of the symmetry currents are given by 

d;A«(x) = (l-^){$(x)$(DW)(x)-$D){x)$(W)(x)} (2-41) 
+ fa Mq {^(x)r 3 r a (D^)(x) - (^Dj 5 )(x)r 3 r a (j^)(x)} , 

&ffl{x) = (l - ^) {${z)${Ity){x) - @D)(x)ty{x)} (2.42) 
+ |a Mq {ij(x)r 3 r a (D%^)(x) - (^D%)(x)r 3 r a ^(x)} . 

Note that these expressions vanish exactly upon summation over x, so that the 
existence of the currents is guaranteed by the lattice Poincare lemma [^] . The 
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symmetry currents themselves will not be needed in the following. Explicit 
expressions for the massless case can be found in |2lJ, for example, but one 
should keep in mind that the integration of eqs. fl2.41 , 2~42^ ) is ambiguous by 
terms with vanishing divergence. 

The Ward identities now take the generic form 

((<5|(u;)S) O) = (S^u)0) , X = A,V, (2.43) 

where O denotes some product of local composite fields. This is an exact 
identity, as the space-time dependent change of variables can be made in the 
functional integral. However, the variation of the composite fields does not 
exactly transform members of a given multiplet among each other^]. Rather 
there are extra terms, which arise due to the non-trivial space-time structure 
of 75. For instance, the axial variation of the pseudo-scalar density is given by 

5 a A (u)P b (y) = \5 ab u;{y)S\y) + |a^(y)rV 7 5 (\u,D]%4>) (y). (2.44) 

We now first assume that u is non-zero only in a single lattice point x ^ y. 
Eq. Q2.44D then reduces to 

6 A (u,)P b (y) = -\au:{x)^{y)T b T a l5 D{y,x) (7^) (x), (2.45) 

where D(y,x) denotes the kernel of the Ginsparg- Wilson Dirac operator. Lo- 
cality of the Ginsparg- Wilson action then implies that the r.h.s. of eq. ( |2.45| ) is 
exponentially small as long as the distance between x and y is large in lattice 
units. We find that this structure is generic and therefore conclude that the 
bare PCAC and PCVC relations hold up to exponentially small corrections, 
provided all operators in the correlation function keep a large distance (in 
lattice units) from the space-time region where the action is varied. 
We now assume that u is constant in a space-time region R, 

fl, if xeR 
w(x) = < (2.46) 
0, otherwise, 



which leads to 



l5 ab S°(y) + |a 5 Y,Hyh5T b T a D(y,x)(i 5 ^)(x), if y E R, 
-la 5 ^2i,(yh 5 T b T a D( y ,x)(i^)(x), if y £ R, 

x&R 



(2.47) 



2 We thank L. Giusti for drawing our attention to this problem. Although this was known 
in the literature (cf. e.g. it has been overlooked by us in an earlier version of this paper. 
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Again, locality of the action implies that the extra terms are exponentially 
small in both cases, provided one has 

minllz — y\\ » a, ifyiR (2.48) 
min||x — y\\ S> a, if y £ R. (2.49) 

xfR 

Furthermore, we note that the above example of the pseudoscalar density is 
generic as the difference to naive continuum considerations always consists of 
terms involving the lattice Dirac operator. Up to exponentially small correc- 
tions, the bare Ward identities ( |2.43| ) can therefore be written in the continuum 
like form 

/ n \ 



(d;A«(x)-2mP a (x)- i ^ a S\x))^X(x t )J = (2.50) 



AiB 

i\xi£R 

« 4 E( (d;V^x)+2^s 3ab P\x))f[ct>tX^i)) = (2-51) 
xeR \ i=l / 

<n)\ a /An) /„ \ \ Arn) 



-I 

i\xi£R 



E (4%, (Ci(-i) • • • Kite) ■ ■ ■ CL(*»)) . 



where we have assumed the conditions ( 2.48 , 2.49 ) to hold with y replaced by 
Xi, for alH = 1, . . . ,n. We have omitted the subscript (M, a) of the correlation 
functions, and the expansion, 

4t("v; w A ) = - ^ (4 r) )^ - ufl (Ta)^ + 0(^ 2 ), (2-52) 

defines the anti-hermitian generators of infinitesimal SU(2) x SU(2) rotations 
in the representation r. 



3 Equivalence between the renormalized theories 

We discuss under which conditions the identity ( |2.35| ) holds in terms of renor- 
malized correlation functions. Without loss one may restrict attention to 
correlation functions of composite fields which keep a physical distance from 
each other. 
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3.1 Renormalized tmQCD 

So far we have dealt with the theory at fixed lattice spacing a, and it is not 
obvious that the theory has a continuum limit. We assume that infrared diver- 
gences of the correlation functions are properly taken care of e.g. by working in 
a finite volume with a suitable choice of boundary conditions. Note that this 
also ensures analyticity of the correlation functions in the mass parameters. 
However, we expect that our conclusions will be valid more generally. 

In perturbation theory, it has been shown that lattice QCD with Ginsparg- 



Wilson quarks is renormalizable [23], and we shall assume that this remains 
true beyond perturbation theory. While this result has been obtained for mass- 
less QCD, we do not expect any additional complication here, as both twisted 
and standard mass terms can be viewed as super-renormalizable interaction 
terms which do not modify the power counting. 

The entire physical information of QCD is contained in the correlation 
functions of gauge invariant composite fields. Having introduced the lattice 
regularized theory, we may completely avoid gauge fixing, and we may also 
avoid the renormalization of the fundamental fields, which only play the role 
of integration variables. We are thus left with the renormalization of the bare 
parameters of the action, and the renormalization of the composite fields which 
enter the correlation functions. The symmetries of the lattice regularized 
theory are the same as in the continuum except for the continuous space-time 
symmetry being replaced by the symmetry group of the hyper-cubic lattice. 
It then follows that renormalized parameters take the form, 

di = Z g9o> m R = z mm, fi K = Zpfiq, (3.1) 

where the renormalization constants are functions, 

Z = Z(g$,an,n q /n,m/fi), (3.2) 

and [i denotes the renormalization scale. 

Even though chiral symmetry is also broken by the mass terms it is cus- 
tomary to renormalize operators such that the chiral multiplet structure of 
the massless theory carries over to the renormalized theory. Renormalized 
operators then are of the form 

teoS = ^ Oil + ct£L* dk '- dk 4l) . (3-3) 

where d k and dy are the mass dimensions of the fields in the multiplets k 
and k'. The structure of this equation follows from the well-known result of 
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renormalization theory which states that composite fields mix under renormal- 
ization with all fields of equal or lower dimension, which transform identically 
under all the symmetries of the regularized theory. We assume here that 
dy < d k , i.e. either there is no mixing with fields of equal mass dimensions, 
or this has already been taken into account by choosing a basis where the 
renormalization matrix is diagonal at the renormalization scale fi. Hence the 
c-coefficients in eq. (|3.3| ) only multiply operators with lower dimensions, and 
this implies that they cannot depend on the renormalization scale jj |24j . While 
the multiplicative renormalization constants Z^ are of the form ([^), the c- 
coefficients may thus be considered functions of the bare parameters, c/q, am 
and a/j, q . Note that the multiplet structure of the bare theory is respected 
by the assignment of a common renormalization constant to all members of a 
multiplet. 

In the following we shall assume that the renormalization of the theory 
and the composite fields works out along these lines. While this is guaran- 
teed in perturbation theory, the non-perturbative renormalization of power 
divergent operators may require an additional effort. In particular, it may 
be necessary to first implement Symanzik's improvement programme to suffi- 
ciently high order in the lattice spacing a before the power divergences can be 
subtracted in an unambiguous way. 



3.2 A special choice of renormalization scheme 



We would like to identify renormalization schemes where the equation ( 2.35[) 



carries over to the renormalized theory, with renormalized fields of the form 
First of all, we notice that the exact PCAC and PCVC relations of the bare 
theory imply that // q and m can be renormalized by the same renormalization 
constant. We may thus choose Z m = implying a multiplicative renormal- 
ization of the bare radial mass M by the same constant. Furthermore, a is 
not renormalized as it is determined by the ratio of the mass parameters. As 
the effect of the chiral rotation is a change of a, we would like to choose all 
multiplicative renormalization constants independently of this angle, i.e. 

Z = Z(glav,M/n). (3.4) 

A simple example is a mass-independent renormalization scheme which is ob- 
tained by renormalizing the theory in the chiral limit p5| . As we have assumed 
an infrared regularization to be in place, such renormalization schemes do ex- 
ist and it is then obvious that eq. (|2.35| ) holds for the renormalized correlation 



12 



functions of multiplicatively renormalizable operators, viz. 

(n<;k (-«')) (teoJ&On) • • • ^)£l(^)) (MRa _ a0 • (3.5) 



.i=i 



The expectation values in the renormalized theory are defined as in the bare 
theory, except that the bare parameters are expressed in terms of the renor- 
malized parameters Mr and g R . Note that eq. ( |3.5| ) is again an exact identity, 
if this is true for ( 2.35| ) and provided the renormalization constants are chosen 



exactly as specified above (i.e. not only up to cutoff effects). 

The case of power divergent operators is slightly more complicated. In 



general, if eq. (3.5) is to be satisfied the c-coefficients cannot be independent 
of a. For example, the renormalized scalar and pseudo-scalar densities are of 
the form 

i(P R ) a = Zp(iP a + 5 a3 a- 3 c P ), (3.6) 
i(5 R )° = Z s {±S° + a- 3 c s ), (3.7) 

with Zg = Zp, and it is well-known that cp vanishes at // q = 0. 

Rather than being independent of a, the additive counterterms satisfy a 
covariance condition. More precisely, assuming that all multiplicative renor- 



malization constants are of the form (3.4), the requirement that eqs. ( 2.35 ) 
and ( |3.5[) hold simultaneously implies, 



R%(-a') c kk>'ic(9l aM, a-a') = c^( ffo 2 , aM, a)R { ^{-a'), (3.8) 

where only the index B is summed, all others being fixed. In the example of 
the pseudo-scalar and scalar densities one finds the equations 

cp(a)\ _ / cos a' sinaA ( cp(a — a')\ 

(3.9) 



eg (a) ) \— sin a' cos a' J \cs{a — a 

where we have only indicated the dependence upon the angle. It is obvious 
that cp + Cg is independent of a and due to the vanishing of cp(0), the solution 
can be parameterised as follows, 

cp(a) = cs(0) sin(a), cs(a) = cs(0) cos(a). (3.10) 

In other words, if the renormalization problem for the scalar density can be 
solved in standard QCD, the solution at all other values of a is given by these 
equations. 
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3.3 Differential equation in a 



The requirements for the renormalization scheme are best characterised by a 
differential equation in a. We start by considering the third flavour component 
of the axial Ward identity (|2.50D , chose the region R to be the space-time 
manifold itself and obtain the exact identity of the bare theory, 



d_ 

da 



nciM = 

Mm \ i= x I (Mja) 

n 



(M,a) 



(3.11) 



where = (T^) 3 . Note in particular that the normalising factor of the 
path integral is independent of a and hence does not generate a disconnected 
contribution. 

Next we recall that analyticity of the renormalized correlation functions 
in the mass parameters m R and /i R is guaranteed by the infrared cutoff. There- 
fore, also their derivative with respect to a must exist, as the relation 



d_ 

da 







d 



Mr 



<9m F 



(3.12) 



follows directly from the definition of the (renormalized) polar mass coordi- 
nates [cf. eq. (2.33|. Differentiating a renormalized n-point function at fixed 
renormalized parameters, applying the chain rule and using eq. ( |3.11| ), one 
obtains the differential equation, 



v. -l>) (teo£i(*i) • • • wfclw 



On) 



1=1 



E T i?B((0R)a^)-(0R)S( 



(M R ,a) 
R >k n A n 



M (x„) 



1=1 



(Mr, a) 



.(3.13) 



Here the differential operator V a is defined through 



_ def 8 



/ 2 ^ 
+ 7T-2- 

a,9R 



a,Af R 



and the coefficients are given by 
d log Z x 



5a 



X = M,g, ki,... ,kn 



M,g 



(3.14) 



(3.15) 
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With a suitable choice of the irrelevant parts of the renormalization countert- 
erms eq. (3.13) holds exactly at finite lattice spacing. The very existence of the 
differentiated correlation function then leads to a constraint for the additive 
counterterms, viz. 



d_ 

da 



M,g 



Jr,r') 
-kk';AC 



- L AB L kk':BC 



(r,r') T {r 
L kk';AB' L BC 



(3.16) 



As one might expect, integrating eq. ( |3.16| ) from a to a — a' reproduces the 
covariance relation ( |3.8| ). Note however, that eq. ( |3.16| ) holds without assum- 
ing eq. (U). 

To relate the renormalized correlation functions defined at two different 
values of the angle a one just has to integrate the differential equation ( 3.13; ). 
It is possible to formally write down the general solution. However, we do not 
find this particularly illuminating and therefore just mention that the general 
relation is much more complicated than eq. (3J3). The question then arises 
under which conditions the relation does take the simple form ( |3.5| ). For this 
to work out, the l.h.s. of eq. ( |3.13| ) should reduce to the partial derivative with 
respect to a, i.e. one requires 



l M M R 







R 



, 2 d 



R a,M R 
(r») 



^>k n A, 



(M R ,a) 



0. 



(3.17) 



Note that such an equation must hold for all renormalized correlation func- 
tions. Therefore, unless Mr and g R are related in a special way, the lx must all 
vanish, i.e. the multiplicative renormalization constants all take the form fl3.4| ). 
We finally remark that eq. ( 3.13j ) expresses the re-parameterisation invari- 



ance of the theory: an infinitesimal change of quark variables of the form (2.27) 
is compensated by a change of a, Mr and g^. Our derivation makes use of the 
axial Ward identity as the change of variables ( 2.27| ) corresponds to a special 
chiral symmetry transformation. However, more general changes of variables 
can be considered using rather general results of renormalization theory |2 



3.4 The role of the Ward identities 



Eq. (3J3) is the relation between renormalized tmQCD and QCD in its sim- 
plest form. Its infinitesimal version is eq. ( |3.13| ) with vanishing coefficients 
lx (3.15). We have argued that it is possible to renormalize the Ginsparg- 
Wilson regulated theory such that eq. ( |3.5| ) holds, and we have worked out 
the conditions on the renormalization constants in this regularization. Based 
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on universality one thus expects that tmQCD in other, not necessarily chirally 
invariant regularizations can again be renormalized such that eq. ( |3.5| ) is sat- 
isfied up to cutoff effects. At least in perturbation theory this can be proved 
rigorously. 

In view of the formulation of tmQCD with Wilson quarks, we would like 
to emphasise the role of the renormalized tmQCD Ward identities, which are 
of the same form as eqs. ( 2.50| , 2.51). As is well-known, the Ward identities en- 



sure that the renormalized composite fields form chiral multiplets, and fix the 
absolute normalization of the symmetry currents. Furthermore, the mass pa- 
rameters are renormalized inversely to the pseudo-scalar and scalar densities, 
and this implies that one may set 

m R = Mucosa, /i R = Mpsina, (3.18) 

with the angle a which remains unrenormalized. 

We are now going to demonstrate that the renormalized tmQCD Ward 
identities also imply that certain linear combinations of the correlation func- 
tions satisfy the Ward identities of standard QCD. In the following we assume 
the continuum limit has been taken and use a shorthand notation for the 
renormalized correlation functions, 

G ^ (^Rti(y)O cxt ) (3.19) 

\ / (Mr, a) 

where ex t is some product of renormalized fields which are localized in the 
exterior of the finite physical space-time region R. Furthermore we assume 
y £ R and denote by G\ v p and Gs the same correlation functions with 
an insertion of c^(Ap)^, ^(Va)^, (-Pr) and (5r)° at the space-time point 
x, integrated over the region R with respect to x. The renormalized Ward 
identities may then be written in the form 

G% - 2M R {cos(a)G P + ism(a)5 3a G s } = -iT£G, (3.20) 
Gy + 2M R sin(a)e 3afc G P = -iTyG, (3.21) 



where the abbreviation 



w^rtrJ^Mw*)^ (*») 



has been used for X = A, V. 

We now consider the two linear combinations 



cG\ + sG^ - 2 Mr {c cos a + s sin a} Gp = -i (cT| + sT$) G, (3.23) 
cGy, - sG l A - 2M R {c sin a - s cos a} Gp = -i (cT^ - sT^) G. (3.24) 
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Multiplying both sides of the equation by the matrices R^ r \a), setting 

R a G^R%(a)((cl )R )^(y)O cxt ) itr x (3.25) 



and also defining 

G^pj = cG^ + sGy, Gyi = cGy — sG]^, (3.26) 

we then find 

RaG\i — 2Mr {c cos a + s sin a} R a Gp = —iR a (cT^ + sTy) G, (3.27) 
RaGy/ - 2M R {c sin a - s cos a} RGp = -iR a (cT^ - s7|) G. (3.28) 

We notice that in order to preserve the canonical normalization of the primed 
currents, one needs c 2 + s 2 = 1, i.e. we may set 



c 



cos(a'), s = sin(a'). (3.29) 



It is then clear that the standard Ward identities can be obtained by choosing 
a' = a, and provided the equations 



^HH^^L + ^^v^cj^k-") = (Tf)^ (3.30) 
flS(«)(c(4 r) )^-s(^)^)41,(-«) = (4 r) )L> (3-31) 



hold. By differentiating with respect to a and using the Lie algebra of the 
SU(2) x SU(2) generators, 

[TlT b A ] = e abc T v , [Tl, I*] = e abc T c A , [T$, T v ] = e abc T$, (3.32) 

one arrives at the conclusion that eqs. (3.3C.3.31) are indeed satisfied. Note 
that the same procedure applies to the remaining flavour components of the 
Ward identities. Hence, the validity of the tmQCD Ward identities implies 
that 

• there exist particular linear combinations of correlation functions which 
satisfy the standard QCD Ward identities for two degenerate quark 
flavours with mass Mr, 

• the linear combinations only depend on the angle a, which is determined 
by the ratio of the quark mass parameters which appear in the Ward 
identities. 
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In particular it is clear that the angle a has no physical significance. We may 
start with any value of a and still obtain the standard chiral flavour Ward 
identities. As a given theory is identified by its symmetries this implies the 
equivalence of theories defined at different values of a, provided the remain- 
ing symmetries are either a-independent or transform covariantly. This is 
certainly true on the level of the renormalized composite fields: for any trans- 
formation of the composite fields at a = one may identify the corresponding 
transformation in the renormalized twisted theory. 

3.5 Concluding remarks 

In practical applications one would like to work with tmQCD at a given fixed 
value of a, and just invoke eq. ( |3.5| ) in order to interpret the results in terms of 
standard QCD. According to the above discussion one may start by imposing 
the tmQCD Ward identities in the renormalized theory. Besides defining the 
value of a this procedure restores the chiral multiplet structure of the bare 
theory. One then still needs to impose a renormalization condition per chiral 
multiplet. If this is done either in the chiral limit or independently of a, 
eq. (|3.5| ) provides the relation to the theory defined at any other angle a, 
including a = 0. The simplification in the Ginsparg- Wilson regularization 
consists in the validity of bare continuum-like tmQCD Ward identities, and 
in the related fact that the bare Ward identity masses coincide with the bare 
mass parameters of the action. Finally we stress that the classical continuum 
theory allows to infer the relation between renormalized theories and 

may hence be used as a guide. 

4 Twisted mass QCD with Wilson quarks 

In this section we discuss in some detail the regularization with Wilson fermions, 
including some practical aspects of applications. 

4.1 Symmetries of the bare theory 

With Wilson quarks the tmQCD Dirac operator is as given in eq. flU]) with 
the usual massless Wilson-Dirac operator^ 

3 

3 For unexplained notation and conventions we refer to ref. |2^] 
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For simplicity we defer the discussion of 0(a) improved tmQCD to a separate 
publication [32]. Here we note that the Wilson term is not left invariant 
by the axial rotation (2.2), and the lattice regulated theories at /i q = and 
/iq 7^ are thus different. This is of course welcome as otherwise the zero mode 
problem would be present in both cases (cf. sect 1). One may think of more 
general lattice Dirac operators, also including a chirally twisted Wilson term. 
However, a moment of thought reveals that this is not really more general, as 
an axial rotation (|2.2| ) may then be used in the lattice theory to eliminate the 
extra Wilson term. Modulo a more general coefficient of the standard Wilson 
term and with re-defined bare mass parameters, one then obtains again the 
action corresponding to eq. (1.1). 

As compared to the theory with JVf = 2 standard Wilson quarks (// q = 0) 
we find that the exact U(2) symmetry is reduced to a U(l) symmetry leading to 
fermion number conservation, and a vectorial U(l) symmetry with generator 
r 3 /2. Concerning the space-time symmetries, the lattice action is invariant 
under axis permutations, whereas reflections such as parity are a symmetry 
only if combined with either a flavour exchange 

F : ^tV, 4>->4>t 1 , (4.2) 

or a sign change of the twisted mass term (i q — > — fi q . We will refer to the 
thus modified parity symmetries as Pp and P respectively. The list of exact 
symmetries is completed by charge conjugation, and we note that twisted mass 
lattice QCD with Wilson quarks has a positive self-adjoint transfer matrix [32|. 



4.2 Renormalized parameters 

As in sect. 3 we assume that infrared divergences are regulated e.g. by working 
in a finite space-time volume with suitable boundary conditions. This implies 
analyticity of the theory in the mass parameters and it is then rather obvi- 
ous that twisted mass lattice QCD is renormalizable by power counting p8| . 
The counterterm structure follows from the symmetries of the regularization. 
Based on these symmetries one concludes that tmQCD is finite after the usual 
renormalization of the coupling and the standard mass parameter, 

gl = Z g gl, m R = Z m m q , m q = m - m c , (4.3) 

and, in addition a multiplicative renormalization of the twisted mass parame- 
ter, 

H K = Z^flq. (4.4) 

In particular we note that the modified parity symmetry, Pp, is sufficient to 
exclude a counterterm oc tr {F^F^}. 
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4.3 Ward identities and renormalization of composite fields 

Following sect. 3 we require that the renormalized theory satisfies the axial and 
vector Ward identities, i.e. the renormalized analogues of eqs. ( 2.50|J2.5T ). Con- 



cerning the vector Ward identity the situation is the same as with Ginsparg- 
Wilson fermions, due to the exact flavour symmetry of massless Wilson quarks. 
Therefore eq. (2.51) holds exactly, with the point-split vector current 



W = ^(x)(^-l)^U{x,^(x + ap) 

+ + aftfru + l)$U{x, ^rV(z)}, (4-5) 

and the local pseudo-scalar density. It then follows that the vector current 
is protected against renormalization, i.e. Zy = 1. More generally, the mul- 
tiplicative renormalization constants of composite fields which belong to the 
same isospin multiplet must be identical in order to respect the vector Ward 
identities. An example is the renormalized pseudo-scalar density which has 
the structure 

i(P R ) a = Z P {iP a + 5 a3 a~ 3 c P } (4.6) 

and cp vanishes exactly at fi„ = 0. The vector Ward identity here implies that 
Zj> is the same for all flavour components, and, moreover, 

Zp = Z-\ (4.7) 

In contrast, the axial Ward identity does not hold in the bare theory. Axial 
Ward identities therefore provide normalization conditions which determine 
finite renormalization constants such as Z\, or finite ratios of scale dependent 
renormalization constants, such as Z$/Zp J|]. Moreover, these finite renor- 
malization constants only depend on the bare coupling <?o and may therefore 
be determined in the chiral limit using standard procedures . Note that the 
finite renormalization constants restore chiral symmetry of the bare theory up 
to cutoff effects. Once this is achieved the renormalization of multiplicatively 
renormalizable fields is similar to the Ginsparg- Wilson case, i.e. the renormal- 
ization constant for a given multiplet is determined by imposing a renormal- 
ization condition on one of its members. Of particular practical interest are 
mass-independent renormalization schemes, which are obtained by imposing a 
renormalization condition at /i q = m q = p5| ]. Based on universality we then 
expect that the relations ( |3.5| ) between renormalized correlation functions hold 
up to cutoff effects. According to sect. 3 the same can be achieved by impos- 
ing a-independent renormalization conditions, where, in the case of Wilson 
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fermions the angle a must be denned through the Ward identity masses (cf. 
subsect. 4.4). 

In principle the Ward identities also determine additive renormalization 
constants which arise due to the explicit breaking of chiral symmetry. An 
example is the renormalization of the iso-singlet scalar density, which has the 
same structure as in eq. ( |3.7D , however, with a coefficient cs which does not 
vanish in the chiral limit. Therefore, the renormalization of the third axial 
Ward identity ( gjg ) requires the explicit subtraction of power divergences. 



While power divergent renormalization problems do not present any particular 
difficulty in perturbation theory, it is less clear how to proceed in a non- 
perturbative approach. A general discussion of this topic is beyond the scope 
of this work. Here we just note that the renormalization of the third axial 
Ward identity may in fact be avoided if one assumes that the physical isospin 
symmetry is restored in the renormalized theory. In the following we will make 
this (plausible) assumption and not discuss the third axial Ward identity any 
further. 

4.4 Definition of the angle a 

According to section 3.4 the angle a is uniquely defined through the renormal- 
ized Ward identity masses. Assuming that Zp has been fixed, the renormalized 
twisted mass parameter is determined due to eq. (^7]) . The renormalized axial 
current and the first two components of the pseudo-scalar density may then 
be used to define m R through the renormalized PCAC relation 

d M (A R )£ = 2m R (P R ) a , a = 1, 2. (4.8) 

In practice one first defines a bare PCAC mass m from the ratio of correlation 
functions involving the bare axial current and pseudo-scalar density. The 
renormalized PCAC mass is then given by 

m R = Z^Zp 1 m = Z m m q . (4.9) 

Using eq. (p~7|), the angle a is then determined as 

tana=^ = -^ = — (4.10) 
m R Zpjn Z m Zpm q 

In general one thus needs the bare PCAC mass m and the axial current nor- 
malization constant Z& to obtain a. Note that the latter is not needed in the 
special case m = 0, i.e. if a = tt/2. 
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4.5 Avoiding lattice renormalization problems 

Once tmQCD has been renormalized in the way described above, eq. ( |3.5| ) can 
be applied to establish a "dictionary" between the renormalized correlation 
functions in QCD and tmQCD. For example, the 2-point function of the axial 
current and the pseudo-scalar density in standard QCD translates as follows 

((A R )l(x)(P R )\y)) {MBfi) = cos(a)<(^ R )J(x)(P R ) 1 (y)> (MR]Q) 

+ sm(a){V 2 (x)(P R ) l (y)) {MRay (4.11) 

More generally, relations between the renormalized composite fields can be 
inferred from the corresponding relations in the classical theory (cf. sect. 2). 
In particular, the above example follows from eqs. ( 2.1C| - 2.12 ). 



As tmQCD and standard QCD with Wilson quarks are not related by a 
lattice symmetry, the counterterm structure for composite fields with the same 
physical interpretation depends upon a. Given a physical amplitude it may 
hence be possible that a particular choice of a leads to simplifications. An 



obvious case is the computation of F n from the 2-point function ( 4.11 ). While 
the standard approach (i.e. the direct computation of the l.h.s.) requires to 
first determine the renormalized axial current, the r.h.s. of this equation at 
a = ir/2 only contains the vector current which is protected against renormal- 
ization. 

Even more interesting is the application of tmQCD to matrix elements 
of the iso-singlet scalar density. At a = tt/2, the physical scalar density is 



represented by the third component of the pseudo-scalar density, see eq. ( 2. 13| ) 



While the scalar density is cubically divergent even in the chiral limit, the 
pseudo-scalar density has a quadratic divergence which vanishes exactly at 
fj, q = 0. The situation is therefore comparable to the case of the renormalized 
scalar density in standard QCD with Ginsparg- Wilson fermions. 

4.6 Inclusion of heavier quarks 

It is straightforward to generalise tmQCD to include any number of heavier 
quark flavours. For the latter one may use the standard regularization with 
(improved) Wilson quarks, as the zero mode problem is practically absent 
at mass parameters which correspond to the physical strange quark mass. 
The renormalization procedure can again be carried out such that the chiral 
flavour Ward identities are respected, which now also involve mixed operators 
of light and heavy quarks. Hence we expect that the "dictionary" between 
tmQCD and standard QCD can again be established by naive continuum 



considerations. As an example we consider the tmQCD continuum action (2.1) 
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and add the action of the strange quark, 

S F [ip,ij,s,s] = J d A x (^+m + ^ q 7 5 r 3 ) ^ + s {p+ m s ) s} . (4.12) 
Using a physical notation, 

^ = ' ^ = ( fi J ) ' ^ 4 ' 13 ) 

the standard PCAC relation, 

9^((?7^75s) = (m' + m s ) e?7 5 s, (4.14) 

is obtained with the rotated axial current and pseudo-scalar density, 

d'^i^s = cos(|a)(i7 A1 75S + zsin(ia)(i7 M s, (4.15) 
^'755 = cos(|a) J75S — i sin(ia) Js, (4-16) 

and with the angle a and the light quark mass m! as given in sect. 2.1. 

4.7 Application to the AS = 2 effective weak Hamiltonian 

Also in the case of operators involving light and strange quarks certain renor- 
malization problems of standard Wilson quarks can be circumvented. An 
interesting example is the AS* = 2 part of the effective weak Hamiltonian, 

AS = 2 = {s llM (l- l5 )d} 2 . (4.17) 

In phenomenology one is mainly interested in the hadronic matrix element of 



this operator between Kq and Kq states [29]. As parity does not change in 
this transition, only the parity conserving part of the operator contributes. 
Hence one decomposes the operator into parity even and odd parts, 

AS = 2 = 0vv+AA - 20 VA . (4.18) 

In the regularization with Ginsparg- Wilson fermions the operator O as=2 and 
thus both Ow+AA and Ova are renormalized multiplicatively. With Wilson 
quarks, the remaining symmetries imply that Cfyv+AA mixes with four other 
parity even operators of the same mass dimension, whereas Ova is still renor- 
malized multiplicatively In tmQCD we now observe that the parity even 
operator in the standard basis is represented by the combination 

°VV+AA = cos(a)Ovv+AA - 2zsin(a)OvA- (4.19) 
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In particular, at a = ir/2, only Ova appears on the r.h.s., and one concludes 
that matrix elements of the physical operator O vv+AA can be computed in 
tmQCD without solving the complicated renormalization problem for the par- 
ity even operator^]. In particular, the Kq — Kq mixing amplitude could be 
extracted from the 3-point function involving Ova and appropriately rotated 
interpolating fields for the kaons [cf. eqs. ( [4.15 , 4~16l )1. 



One might be worried that additional counterterms to Ova may be nec- 
essary in tmQCD. As there is no such term at [i q = 0, possible counterterms 
must be accompanied by at least one power of the twisted mass parameter, 
and the flavour structure requires them to be again four-quark operators. For 
dimensional reasons such counterterms can only contribute cutoff effects of the 
order afi q , and a closer look shows that the parity even operators multiplied 
by a/i q are indeed allowed by the tmQCD symmetries. 

4.8 Technical complications 

The equivalence between tmQCD and standard QCD is a statement about the 
renormalized theories in the continuum limit. When tmQCD is used to define 
the standard QCD correlation functions, some of the physical symmetries are 
only restored in the continuum limit. In particular, this is the case of the 
flavour symmetry and parity, which are exact lattice symmetries in standard 
QCD with Wilson fermions, but which are only recovered in the continuum 
limit if a ^ 0. In practice the problem shows up e.g. as an ambiguity in the 
definition of a, which is induced by the usual ambiguity in the definition of 
the critical mass m c by terms of 0(a) (or 0(a 2 ) if the theory is improved) p7fl . 
To illustrate the consequences consider the computation of the pion mass at 
fixed cutoff a, using the 2-point function 

G{xq - yo) = a 3 ^2(P'\x)P' 3 (y)). (4.20) 



In tmQCD this correlation function is computed by replacing P' as in eq. ( 2.12j ), 



where the relative multiplicative renormalization of S° and P 3 is assumed to 
have been fixed by the axial Ward identity. Note also that the exponential 
decay of this correlation can be determined without knowledge of the additive 
renormalization constants. 

It is now obvious that both the ambiguity in a and the 0(a) ambiguity in 
the relative renormalization of the densities imply that the correlation func- 
tion ( |4.20| ) contains cutoff effects which are proportional to the propagator of 
the (physical) scalar density S'°. An analysis of the 2-point function (4.20) 



4 



For an alternative proposal see ref. |31 
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at fixed a must therefore take into account the states with the quantum num- 
bers of an iso-singlet scalar state. We note in passing that 0(a) improvement 
does not alter this situation, as it merely reduces the ambiguities to 0(a 2 ). 
In particular, at fixed cutoff, the relevant symmetries remain inexact and the 



qualitative behaviour of the 2-point function ( 4.2C ) remains unchanged. 

In general the analysis of the hadron spectrum in tmQCD at fixed a 
must include states which are allowed by the lattice symmetries but have the 
"wrong" continuum quantum numbers. Although this is not a fundamental 
problem, the analysis is somewhat more complicated than in lattice QCD with 
standard Wilson quarks. We also note a side effect for the determination of 
hadronic matrix elements. When these are extracted from tmQCD correlation 
functions it may not be necessary to increase the distances until the desired 
physical state is completely isolated. It is sufficient to establish that contri- 
butions from the excited states with the correct continuum quantum numbers 
are negligible, as all other contaminations merely modify the cutoff effects of 
the matrix element. 

5 Conclusions 

In this paper we have advocated the use of twisted mass QCD with Wil- 
son quarks as an alternative regularization of QCD with two degenerate light 
quarks. Using Ginsparg- Wilson fermions as a tool, we have demonstrated in 
what sense tmQCD is equivalent to standard QCD. In particular, we have 
clarified under which conditions the relations between renormalized correla- 
tion functions take the simple form fl3.5| ), which is the quantum analogue of 
the naive relations derived in the classical continuum theory. 

Twisted mass lattice QCD provides a clean field theoretical solution to the 
problem of unphysical zero modes. While our work on tmQCD is motivated 
by this problem, we also observe a few additional benefits. In particular we 
have given examples where renormalization problems of lattice operators can 
be circumvented by working in the fully twisted theory with a = tt/2. For the 
sake of simplicity we did not discuss 0(a) improvement of tmQCD. This topic 
is deferred to a separate publication ]32| . 

First numerical simulations using (quenched) tmQCD have already been 
carried out, and a scaling test in a small volume has been presented in |p3| . 
It is hoped that the chiral limit can be approached much more closely in 
tmQCD than previously possible with Wilson quarks. In particular the AL- 
PHA collaboration plans to extend the work of [34] to much smaller quark 
masses where (quenched) chiral perturbation theory should be safely applica- 



ble [35]. Furthermore, a project to determine the Kq — Kq mixing amplitude 
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using tmQCD is underway [36|. In the future, it will also be interesting to see 
whether numerical simulations of full QCD can benefit from using a twisted 
mass term. 

This work is part of the ALPHA collaboration research programme. We 
would like to thank M. Delia Morte, M. Liischer, G.C. Rossi, R. Sommer and 
A. Vladikas for useful discussions. S. Sint acknowledges partial support by 
the European Commission under grant no. FMBICT972442. P.A. Grassi is 
supported by the National Science Foundation (grants no. PHY-9722083 and 
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